Condensation energy in Eliashberg theory — from weak to strong coupling 
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We consider two issues related to the condensation energy in superconductors described by the 
Eliashberg theory for various forms of the pairing interaction, associated either with phonon or 
electronic mechanisms of superconductivity. First, we derive a leading correction to the BCS formula 
for the condensation energy to first order in the coupling A. Second, we show that at a given A, the 
value of the condensation energy strongly depends on the functional form of the effective pairing 
interaction F(oj). 
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A non-monotonic doping dependence of superconduc- 
tivity in high T c cuprates revived the interest to the issue 
of the condensation energy in strongly coupled supercon- 
ductors P, 0, IE S IE @ ■ The condensation energy E c is 
the energy gain in a superconducting state compared to 
a normal state at the same T. In a BCS superconductor, 
E c smoothly increases below T c and at T = reaches 
e bcs _ _t/tv A 2 /2, where V is the volume, A is the 
superconducting gap and N is the fermionic density of 
states 0- The decrease in the total energy upon pair- 
ing results from a fine competition between an increased 
kinetic energy and a decreased 

potential energy, both of which are much larger than 
E c . Within BCS theory, the condensation energy is re- 
lated to the jump of the specific heat at T c as C s — C n ~ 
%mE c /T c . 

The BCS formula for the condensation energy is valid 
at weak coupling. Experimental discoveries of strong 
fermionic self-energy in the normal state of the cuprates 
stimulated the studies of the condensation energy away 
from the BCS limit. In earlier papers |fj Haslinger and 
one of us analyzed the condensation energy vs coupling 
for spin-mediated d— wave superconductivity. It was 
found numerically that the condensation energy rapidly 
deviates from the BCS form as coupling increases, then 
saturates and then even decreases at strong coupling, due 
to feedback on the bosonic spectrum from the pairing. 
The crossover from weak to moderate coupling, however, 
has not been analyzed in detail in Ref 6] . 

In the present paper, we consider condensation en- 
ergy in superconductors described by the Eliashberg the- 
ory, for various forms of the pairing interaction T(lu) oc 
(wjj+w 2 ) -7 ' 2 (0 < 7 < 2), associated either with phonon 
or electronic mechanisms of superconductivity (we dis- 
cuss the examples below). We address two issues that 
were not discussed in Ref. jg. First, we analyze, both 
analytically and numerically, the leading correction to 
the BCS form of the condensation energy at weak cou- 
pling. We show that the relative correction scales as a 
dimensional coupling A and does not depend on the func- 
tional form of the pairing interaction (i.e., on 7). Second, 
we compare the condensation energy for fixed A, as a 
function of 7. In the limit of A — > 00, E c is the largest 
at 7 = 2, and for this 7 it actually diverges as log A. 



However, we found that this behavior only holds for ex- 
tremely large A ~ 10 2 — 10 3 . At moderate A, we found 
numerically an opposite trend: a strong increase of E c 
with decreasing 7. We understood this behavior analyt- 
ically as originating already within the BCS theory and 
coming from the strong 7 dependence of the prefactor of 
the superconducting gap. 

We consider a model in which low-energy fermions are 
interacting by exchanging either phonons or collective 
spin or charge fluctuations. This exchange gives rise 
to an attraction in one of pairing channels, such that 
at T = 0, the system is in the superconducting state. 
At weak coupling, both phonon-mediated and collective 
mode-mediated pairings can be treated within BCS the- 
ory. At moderate and strongcoupling, we assume that 
the Eliashberg theory is valid Q, i.e., that the self-energy 
is large, but it predominantly depends on frequency, and 
vertex corrections can be still neglected. Eliashberg the- 
ory can be straightforwardly justified for phonon super- 
conductors due to the smallness of the the sound velocity 
compared to the Fermi velocity 0, For electronic 

mechanism of pairing, it is justified on different reasons, 
but still, at low energies bosonic excitations turn out to 
be slow modes compared to electrons, i.e., the effective 
bosonic velocity is smaller than vp. We discuss the jus- 
tification of the Eliashberg theory for electronic pairing 
mechanism in some detail at the end of the paper, after 
we discuss the results. 

The smallness of (real of effective) bosonic velocity 
compared to vp implies that in the pairing problem, 
the momentum integration is factorized: the one over 
momenta transverse to the Fermi surface involves only 
fermions, while the one along the Fermi surface involves 
the bosonic propagator x(Q\\ ^Qperp = 0,u>). As a con- 
sequence, the theory operates with the effective "local" 
pairing interaction T w cx f dqux(q\\,oj). This interaction 
quite generally can be cast in the form [TT| 



(i+w7<4) 



- 7 /2 



(1) 



where LOp> is characteristic frequency. The case 7 = 2 cor- 
responds to phonon-mediated superconductivity, where 
ojp> is a Debye frequency. Other values of 7 correspond 
to electronic pairing. For instance, 7 = 1/3 corresponds 



2 



to the pairing mediated by 2D ferromagnetic Ising spin 
fluctuations or long- wavelength charge fluctuations |12| . 
In both cases, x(qS w ) * (£~ 2 + q 2 + M/ 1 ?)- Integrat- 
ing over q, one immediately finds T x — > const at x — > 
and T x oc (x/£ _3 ) -1 / 3 at large x, in agreement with |JTJ 
for 7 = 1/3. The behavior at intermediate x is more 
complex than in (JIJ but this does not affect the physics. 
Similarly, 7 = 1/2 corresponds to the pairing mediated 
by 2D overdamped spin or charge fluctuations peaked at 
a finite momentum |l3l , while 7 = 1 corresponds to 
the pairing mediated by propagating spin or charge den- 
sity waves. The case 7 — ► describes all 3D pairings 
mediated by either charge or spin fluctuations . 

The variables of the Eliashbcrg theory are the pair- 
ing gap A(w) and the quasiparticle renormalization fac- 
tor Z(u>). They are related to the fermionic self-energy 
£(fc,w) and the pairing vertex $(fc,w) by uj + Y,(uj) = 
u>Z(u>), A(u) = $(u)/Z(u)). The equation for A is 
decoupled from that for Z for any form of T w and 
is obtained by straightforward extension of that for 
phonons ^3 
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The dimensionless A is defined as A = (g/ojn) 1 where g 
is the effective fermion-boson interaction. This definition 
implies that at ujd — > 0, the gap equation becomes in- 
dependent on u>d as it indeed should be as ujd — * just 
implies that a pairing boson becomes gapless, i.e., the 
pairing occurs near a quantum critical point (QCP). As 
in previous studies |lf)j | , we fix g to set the overall energy 
scale, and compute E c as a function of A. The choice to 
measure E c in units of g is justified on the grounds that 
near QCP, g does not critically depend on the distance 
to criticality. 

Once A w , is found by solving |J2j, Z^ , is immediately 
obtained from 

00 

Z Ujn , = l + nTX J2 1 a ' — r Wm ,-^ (3) 

The expression for the condensation energy in the 
Eliashberg theory for phonons (7 = 2) has been obtained 
by Wada and Bardeen and Stephen Rja. The exten- 
sion to arbitrary 7 is straightforward [ll|. In terms of 
Du> m = A ^ m and T Wm the condensation energy reads 



E r . 
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If we approximate A w in the first term by a constant, this 
term yields a BCS condensation energy which at T = 
is E Ct Bcs — —NqA 2 /2. The combination of the second 
and the third terms in (J3J , and the frequency dependence 
of Am, account for the corrections to BCS formula. We 



analyzed both corrections and found that at small A, the 
leading correction comes from the combination of second 
and third terms. At T = 0, the sum of these two terms 
reduces to 



5E C 



N X 







(1 + ((w - W ')A^) 2 ) 7/2 (1 + (( W + w')M>) 2 ) 7/2 , 



(5) 



Simple considerations show that the dominant contribu- 
tion to the double integral comes from the range when 
either u ~ u>£>, and a/ ~ A, or vice versa. Since at weak 
coupling, A << cJd, one can expand in the smaller fre- 
quency in the last term in j£| . One then obtains from 



7 



dujuj 



5E c = -N X^ r , 

"I Jo (1 + {u/ujdY) 1 ^' 2 



aui OJ — ; 



(6) 
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The last integral is the same one that yields a BCS con- 
densation energy. Evaluating the first integral, we obtain 



such that 



SE C = XE c bcs 



E c = -No— (1 + A) 



(J) 



(8) 



Note that the relative correction to the BCS result de- 
pends only on A, but not on 7. 

Eq. © for E c can be formally obtained if one assumes 
that the quasiparticle renormalization factor Z(u>) can 
be approximated by a constant, and this constant does 
not change between the superconducting and the nor- 
mal state, where Z = 1 + A 10]. For a constant Z, the 
full Green's function differs from the BCS result only by 
the renormalization — > e* k — e^/Z. Since the con- 
densation energy is the integral over of the product 
of the Green's function and the self-energy, the effect of 
Z can be absorbed into dek = Z del. Obviously then 
E c = E Ci bcs (1 + A) as in Eq. JSJ 10]. It is not a'priori 
clear, however, whether the change in Z between the nor- 
mal and the superconducting state can be neglected, as, 
by rough estimates, the contribution to E c from Z sc — Z n 
is of the same order as in Q . We checked this explicitly 
for 7 = 2 using the Bardeen-Stephen formula for E c and 
found that the correction to E c due to Z ac — Z n is small, 
because of an extra frequency integration involved, and 
scales as XE c ^bcs * (A/Xud) 2 "C XE c _bcs- Our explicit 
calculation of SE C above shows that this smallness sur- 
vives for all 7. 

To test this result, we computed E c numerically for a 
range of small to intermediate A. 

The equation for A UJm was solved using an iterative 
computational method. Two representative sets of re- 
sults for 

Ao; m , Z Um , So, ro and for A = 0.3 and 7 = 2 and 
7 = 1 are shown in Fig^ The condensation energy vs 
A is plotted in Fig [2 At the lowest A < 0.5, we indeed 

A 2 

reproduced the BCS result E Ct scs = —No — . We see, 
however, that at intermediate A < 1, the condensation 
energy deviates from the BCS formula, but closely follows 
Eq. (jHJ. At even larger A, E c deviates even from the 
modified BCS expression. 

We next discuss the behavior of the condensation en- 
ergy at large A. Here the dependence on 7 becomes cru- 
cial. First, we found from numerical computation that for 
7 = 2, the condensation energy never saturates and keeps 
slowly increasing with A (see Fig|3]and Fig^J. This can 
be easily understood analytically: using lud — g/vX and 
substituting Ar Wm = + A _1 g 2 ) into Eq. Q, one 

obtains that the second term in (JIJ diverges as log A . 
A more careful examination shows |ll| that this diver- 
gence actually comes from the divergence of the free en- 
ergy of the normal state [l^ 



F n = N Q ir 2 T 2 g 2 



E 



1 - sgn(u) m )sgn(u r , 

(W m - UJm') 2 + X' 



lg2 



(9) 
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FIG. 1: The results for (a) E Wm) (b) Z Um , (c) and 
(d) A ulm . Upper panel, 7 = 2, lower panel 7 = 1. In both 
cases, we set A = 0.3. The dashed lines in the plots for E Um 
indicates that E„ m is linear in frequency at small ui m . The 
insert in the figure for A„ m shows that the gap is nearly a 
constant for small u m . 
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FIG. 2: The condensation energy, E c . Dashed, dotted and 
solid lines represent the BCS result, Eq. (6), and our numer- 
ical result for E c respectively. 



At T = and A — > 00 (i.e., ujd — * 0), the 2D integral 
in © diverges as log A. 

At the same time, in the superconducting state, the 
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FIG. 3: The condensation energy — E c for 7 = 2 com- 
puted numerically for large A (solid line), and the analytical 
Ec n (\) <x log A within the same range (dashed line). 



free energy 

F sc = 2ttTN ^ 



ir 2 T 2 N \ x 



A2 



W m W m ' + Aa> m A a 



A 2 



is free from divergences as one can easily verify. Obvi- 
ously then, E c = F sc — F n diverges as log A. In FigElwc 
plotted E c at large A vs log A. 

We clearly see logarithmic behavior at A ^> 1. 

The results for other 7 are presented in Fig0] For 
7 < 2, the condensation energy does not diverge at 
A = 00 as the frequency integral in F n in Eq. J3J now 
converges. This implies that at large A, E c decreases by 
magnitude when 7 decreases from 7 = 2. We indeed 
reproduced this behavior in our numerical calculations 
(see the insert in Fig0J. We, however, found numerically 
that E c (j = 2) becomes the largest by magnitude only 
at enormously large A ~ 10 3 . At moderate couplings, we 
actually found an opposite trend - \E C \ increases with de- 
creasing 7, and the increase becomes quite strong when 
7 becomes smaller than 1 (see the bulk of Fig. 0J. This, 
e.g., implies that for same A > 1, in systems where the 
pairing is mediated by propagating optical phonons, the 
condensation energy in 3D is larger than in 2D. It also 
implies that, at the same moderate coupling A, the 2D 
pairing by ferromagnetic (Ising) spin fluctuations should 
yield larger E c than 2D pairing by Heisenberg antifer- 
romagnetic spin fluctuations, and that for the later, the 
condensation energy at a given A is larger for Ornstein- 
Zernike form of the spin susceptibility near (tt, tt) than for 
flat spin susceptibility near (71% tt) (this flat susceptibility 
emerges in RPA-based studies based on the measured 
Fermi surface in Bi2212 |!s|). 

To understand the origin of this effect, we recall that 
for moderate couplings, one can use the modified BCS 
formula, Eq. JSJ. Since 7 enters this formula only 
through A, large difference in E c between different 7 
must be due to a difference in the values of A. To esti- 
mate the 7 dependence of A at a given A ~ 1, we just 
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FIG. 4: The condensation energy for various 7. The insert 
shows the large- A behavior of E c . At weak coupling, \E C \ 
increases with decreasing 7, while at very strong coupling 
it becomes the largest for 7 = 2. Observe that E c reaches 
saturation faster as 7 decreases. 



solve the gap equation to a logarithmic accuracy. Ap- 
proximating A w by a constant Aq, we obtain from Eq. 



(2) 



1 
A 



dx- 



1 



x(x 2 + !)'■ 



(11) 



where Ao = Aq/ljd. At 7 — > (i.e., when the pairing 
interaction tends to a constant), the integral diverges at 
the upper limit. This divergence is artificial as in real- 
ity the integration must not go beyond x max ~ Ep/^D- 
We, however, will focus on finite 7, when the integral 
converges. Carrying out the integration, we obtain an 
effective BCS formula 



A 



uj d e 



exp 



*(1) 



(12) 



where ^(x) is the digamma function. The first term in 
the r.h.s. of l|12fl is the conventional weak-coupling re- 
sult for a phonon superconductor. The dependence on 7 
emerges through the second factor. We found that the 
dependence on 7 is actually quite strong already at in- 
termediate 7. We plot our numerical A(7) obtained for 
A = 0.5 against EalT2l in Fig|3J We see that the agree- 
ment between numerical and analytic Aq is quite good. 
As, to a very good accuracy, the condensation energy for 
A = 0.5 scales as E c oc Aq, the good agreement in Fig. 
|J5J implies that the origin of a very strong variation of E c 
with 7 is indeed a strong 7— dependence of the magnitude 
of A. 

Finally, we discuss the validity of the Eliashberg the- 
ory. For phonon superconductors, Eliashberg theory is 
based on the smallness of the ratio of the electron mass 
and the ionic mass, or, alternatively, with the smallness 
of the sound velocity v s compared to the Fermi velocity 
vf 0- The theory is valid as long as Xv s /vp < 1, where 
A is the dimensionless coupling constant. On the scale of 
this parameter, one can neglect vertex corrections which 
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FIG. 5: A 2 (7) from Eg.lTHl (line), together with the numer- 
ical results for Ao 2 for seven different 7 (dots). We used 
A = 0.3. The inset shows a close view of A 2 (7) for 1 < 7 < 2. 

are small in Xv s /vf- Due to the smallness of v s /vf, the 
dimensionless A doesn't need to be small, and one can 
analyze both weak coupling A < 1 and strong coupling 
1 < A < vf/v s limits. 

For pairing mediated by electronic collective modes, 
Eliashberg theory is justified on different grounds. The 
velocities of collective modes are of the same order as 
the Fermi velocity, hence a direct argument based on 
v s /vp ratio would not work. It turns out, however, 
that when the velocities of the two modes are compa- 
rable, the renormalization of the bosonic propagator due 
to Landau damping mechanism (which was unimportant 
for phonons) becomes relevant, and for A > 1, collective 
modes become diffusive excitations at typical momenta 
and frequencies relevant to pairing. This ensures that 
at small frequencies, bosons become slow modes com- 
pared to fermions. This is what Migdal theorem requires 
on physical grounds, as vertex corrections involve vir- 
tual processes in which fermions are forced to oscillate 
at bosonic frequencies. Once bosons are slow mode com- 
pared to electrons, fermions oscillate far from their own 
mass shell, and vertex corrections are reduced. The cal- 
culations indeed show that at A > 1, vertex corrections 
are much smaller than dY,(uj)/duj |L%ll4| . For example, 
for the pairing by 2D ferromagnetic and incommensu- 



rate spin fluctuations, dY,{uj)/duj = A at large A, while 
vertex corrections remain 0(1) |12|. For the pairing 
by 2D antiferromagnetic spin fluctuations, vertex cor- 
rections grow with A, but only as log A 14] . In both 
cases, there is a wide range of A > 1 where frequency- 
dependent self-energy is already large, i.e., the system 
falls into the strong coupling regime, but vertex cor- 
rections still can be safely neglected. In this regime, 
Eliashberg theory is valid. Note, however, that unlike 
the case of electron-phonon interaction, Eliashberg the- 
ory for electron-electron interaction is only applicable at 
A > 1. 

There are two other complications for collective mode- 
mediate pairing. First, the interaction is generally 
momentum-dependent, and the pairing is not of s — type. 
This does not invalidate Eliashberg theory as vertex cor- 
rections still remain small, and the momentum integra- 
tion in the gap equation is still factorized. As a result, 
Afc i(i) ss AfcA^, and the momentum dependence of the 
gap only affects the functional form of T u , but does not 
change the structure of the equation for A u . Second com- 
plication is that the propagator of a collective mode by 
itself changes once the system enters a superconducting 
state. This, roughly, implies that 7 changes with T/T c . 
In the analysis above we focused on T = 0, assuming that 
all feedback effects are already incorporated into T^. 

To conclude, in this paper we computed the condensa- 
tion energy E c for Eliashberg superconductors for a wide 
range of coupling strengths and for several different phys- 
ically motivated forms of the pairing interaction (param- 
eterized by 7), representing both phonon (7 = 2) and 
electronic collective mode mediated superconductivity 
(0 < 7 < 2). We found analytically and confirmed and 
numerically that for weak coupling, E c — E c ^bcs(^ + A), 
where E c , B cs = -N A 2 /2 is BCS result. ' This result 
holds for all 7. We also found that in the opposite limit 
of very large A, E c diverges as log A if the pairing is medi- 
ated by phonons, but remains finite for electronic mech- 
anisms of pairing and decreases with decreasing 7. How- 
ever, for moderate A > 1, we found an opposite trend - 
the condensation energy at a given A strongly increases 
with decreasing 7. We argued that this behavior holds in 
the regime where E c oc A 2 , and is a consequence of the 
strong 7 dependence of the pairing gap A. We computed 
E c ("f) analytically and numerically and obtained a very 
good agreement between the two results. 

We acknowledge useful discussions with A. Abanov, 
B. Altshuler and E. Yuzbashyan. The research was sup- 
ported by NSF DMR 0240238. 
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